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We present a comprehensive analysis of quantum fluctuation effects in the superfluid ground state of an 
attractively interacting Fermi system, employing the attractive Hubbard model as a prototype. The superfluid 
order parameter, and fluctuations thereof, are implemented by a bosonic Hubbard-Stratonovich field, which 
splits into two components corresponding to longitudinal and transverse (Goldstone) fluctuations. Physical 
properties of the system are computed from a set of approximate flow equations obtained by truncating the 
exact functional renormalization group flow of the coupled fermion-boson action. The equations capture the 
influence of fluctuations on non-universal quantities such as the fermionic gap, as well as the universal infrared 
asymptotics present in every fermionic superfluid. We solve the flow equations numerically in two dimensions 
and compute the asymptotic behavior analytically in two and three dimensions. The fermionic gap A is reduced 
significantly compared to the mean-field gap, and the bosonic order parameter a, which is equivalent to A in 
mean-field theory, is suppressed to values below A by fluctuations. The fermion-boson vertex is only slightly 
renormalized. In the infrared regime, transverse order parameter fluctuations associated with the Goldstone 
mode lead to a strong renormalization of longitudinal fluctuations: the longitudinal mass and the bosonic self- 
interaction vanish linearly as a function of the scale in two dimensions, and logarithmically in three dimensions, 
in agreement with the exact behavior of an interacting Bose gas. 

PACS numbers: 05.10.Cc, 67.25.D-, 71.10.Fd 



I. INTRODUCTION 

Interacting Fermi systems exhibit very diverse behavior on 
different energy scales. Composite objects and collective phe- 
nomena emerge at scales far below the bare energy scales of 
the microscopic degrees of freedom. This diversity of scales is 
a major obstacle to a numerical solution of microscopic mod- 
els, since the most interesting phenomena appear only at low 
temperatures and in systems with very large size. It is also 
hard to deal with by conventional many-body methods, if one 
tries to treat all scales at once and within the same approxima- 
tion, such as a partial resummation of Feynman diagrams. 

Therefore, it seems natural to integrate degrees of free- 
dom (bosonic and/or fermionic fields) with different energy 
scales successively, descending step by step from the high- 
est scale present in the microscopic system. This gener- 
ates a one-parameter family of effective actions which inter- 
polates smoothly between the bare action of the system, as 
given by the microscopic Hamiltonian, and the final effective 
action from which all physical properties can be extracted. 
The vertex functions corresponding to the effective action at 
scale A obey an exact hierarchy of differential flow equations, 
frequently referred to as "exact" or "functional" renormal- 
ization group (fRG).' This hierarchy of flow equations is a 
transparent source of truncation schemes for a given physical 
problem.^ 

Most interacting Fermi systems undergo a phase transition 
associated with spontaneous symmetry breaking at low tem- 
perature. In a renormalization group flow, the instability of 
the normal state is signalled by a divergence of the effective 
two-particle interaction at a finite scale A^. , that is, before all 
degrees of freedom have been integrated out. To continue the 
flow below the scale Ac, the order parameter corresponding 
to the broken symmetry has to be implemented. This can be 



accomplished in a purely fermionic theory or, alternatively, 
by introducing a bosonic order parameter field via a Hubbard- 
Stratonovich transformation. 

In a fermionic description the order parameter can be imple- 
mented by adding an infinitesimal external symmetry break- 
ing field to the bare Hamiltonian, which is then promoted to 
a finite order parameter below the scale To capture first 
order transitions, one has to implement the order parameter 
via a finite counterterm, which is added to the bare Hamilto- 
nian and then subtracted again in the course of the flowii A 
relatively simple one-loop truncation^ of the exact fermionic 
flow equation hierarchy yields an exact solution of mean-field 
models.^''*''' The same truncation with momentum-dependent 
vertices provides a promising approximation to treat sponta- 
neous symmetry breaking also beyond mean-field theory. Re- 
cently it has been shown that this approximation yields quite 
accurate results for the fermionic gap in the ground state of 
the two-dimensional attractive Hubbard model for weak and 
moderate interactionsi^ 

If the order parameter is implemented via a Hubbard- 
Stratonovich field, one has to deal with a coupled theory of 
bosons and fermions. The evolution of the effective action of 
this theory is again given by an exact hierarchy of flow equa- 
tions for the vertex functions. A truncation of this hierar- 
chy has been applied a few years ago to the antiferromagnetic 
state of the two-dimensional repulsive Hubbard model.** Im- 
portant features of the quantum antiferromagnet at low tem- 
peratures were captured by the flow. More recently, various 
aspects of superfluidity in attractively interacting Fermi sys- 
tems have been studied in the fRG framework with a Hubbard- 
Stratonovich field for the superfluid order parameter. Ap- 
proximate flow equations were discussed for the superfluid 
ground state," for the Kosterlitz-Thouless transition in two- 
dimensional superfluids'^ and for the BCS-BEC crossover in 
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three-dimensional cold atomic Fermi gasesji^ 



II. BARE ACTION 



In this work we analyze the fRG flow of the attractive Hub- 
bard model, using the Hubbard-Stratonovich route to symme- 
try breaking. The attractive Hubbard model is a prototype of 
a Fermi system with a superfluid low temperature phase. '"^ 
In particular, it is a popular model for the crossover from 
BCS-type superfluidity at weak coupling to Bose-Einstein 
condensation of strongly bound pairs at strong coupling .i^^ii^ 
An experimental realization of the attractive Hubbard model 
is conceivable by trapping cold fermionic atoms in an op- 
tical lattice and tuning the interaction close to a Feshbach 



resonance.' 
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We focus on the superfluid ground state. The importance 
of quantum fluctuations in the superfluid ground state has 
been emphasized recently in the context of the BCS-BEC 
crossover. Although the long-range order is not destroyed 
by fluctuations in dimensions d > I, the order parameter cor- 
relations are nevertheless non-trivial in d < 3. The Gold- 
stone mode leads to severe infrared divergences in perturba- 
tion theory. A detailed analysis of the infrared behavior of 
fermionic superfluids has appeared earlier in the mathemati- 
cal literature,^' where the perturbative renormalizability of the 
singularities associated with the Goldstone mode was estab- 
lished rigorously. To a large extent divergences of Feynman 
diagrams cancel due to Ward identities, while the remaining 
singularities require a renormalization group treatment. Since 
the fermions are gapped at low energy scales, the infrared be- 
havior of the collective, bosonic sector in fermionic superflu- 
ids is equivalent to the one of an interacting Bose gas, where 
the Goldstone mode of the condensed state strongly affects the 
longitudinal correlations, leading to drastic deviations from 
mean-field theory in dimensions d < 3.22^ 

The purpose of our work is to construct a relatively simple 
truncation of the exact fRG flow which is able to describe the 
correct infrared asymptotic behavior, and which yields rea- 
sonable estimates for the order parameter at least for weak 
and moderate interaction strength. From a numerical solution 
of the flow equations, which we perform in two dimensions, 
we obtain information on the importance of Goldstone modes 
and other fluctuation effects. 

In Sec. II we introduce the bare fermion-boson action ob- 
tained from the attractive Hubbard model by a Hubbard- 
Stratonovich transformation. Neglecting bosonic fluctuations, 
one recovers the standard mean-field theory for fermionic su- 
perfluids, as recapitulated in Sec. III. By truncating the ex- 
act fRG hierarchy, we derive approximate flow equations in- 
volving fermionic and bosonic fluctuations in Sec. IV. At the 
end of that section we reconsider mean-field theory from a 
flow equation perspective. Sec. V is dedicated to a discus- 
sion of results obtained by solving the flow equations. We 
discuss the asymptotic behavior in the infrared limit in two 
and three dimensions and then present numerical results for 
the flow in two dimensions, where fluctuation effects are most 
pronounced. Finally, we summarize our results in Sec. VI. 



As a prototype model for the formation of a superfluid 
ground state in an interacting Fermi system we consider the 
attractive Hubbard model 
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where c. and c. are creation and annihilation operators for 
spin-^ fermions with spin orientation cr on a lattice site i. For 
the hopping matrix we employ t\j - -f if i and j are near- 
est neighbors on the lattice, and t\j - otherwise. On a d- 
dimensional simple cubic lattice, this leads to a dispersion re- 
lation = -2f 2f=i cos A;,-. For the attractive Hubbard model 
the coupling constant U is negative. 

The attractive Hubbard model has a superfluid ground state 
for any particle density n in li > 2 dimensions, provided the 
lattice is not completely filled (n - 2) or empty (« = 0). At 
half filling (« = 1) the usual U(l) global gauge symmetry be- 
comes a subgroup of a larger SO(3) symmetry group, and the 
order parameter for superfluidity mixes with charge density 
order»ii 

Our analysis is based on a functional integral representa- 
tion of the effective action, that is, the generating functional 
of one-particle irreducible correlation functions. For the Hub- 
bard model, the starting point is a functional integral over 
fermionic fields iff and tfr with the bare action 



,0-] = - r ' 

Jkcr 
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+ I U ^-k+i[^k+i,^>Pk'+%^>P-k'+U , (2) 

Jk,k',q ... 

where = - is the single-particle energy relative to the 
chemical potential. The variables k = (ko, k) and q - (qo, q) 
collect Matsubara energies and momenta. We use the short- 
hand notation = ^ = ^ for momentum 

and energy integrals, and J^^ includes also a spin sum. We 
consider only ground state properties, such that the energy 
variables are continuous. 

The attractive interaction drives spin-singlet pairing with 
s-wave symmetry and a spontaneous breaking of the global 
U(l) gauge symmetry. Therefore, we decouple the Hubbard 
interaction in the s-wave spin-singlet pairing channel by in- 
troducing a complex bosonic Hubbard-Stratonovich field (pg 
conjugate to the bilinear composite of fermionic fields^! 



(3) 



This yields a functional integral over iff, t// and <p with the new 
bare action 

Foliff, 1^, 0] = - (ffkAiko - ^k) 4fkiT - J^'^*/'^'^'? 

+ I {'k-k+'il4'k+i\^q + 4'k+'it<P-k+'ilfq)- (4-) 
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where ^* is the complex conjugate of 0, while t/r and i/r are 
algebraically independent Grassmann variables. 

Our aim is to compute fermionic and bosonic correlation 
functions with a focus on a correct description of the low- 
energy (infrared) behavior. A central object in our analysis is 
the effective action r[t//, 0], which can be obtained by func- 
tional integration of the bare action in the presence of source 
fields and a subsequent Legendre transform with respect to 
these fields. Functional derivatives of r[i/', i^, 0] with respect 
to i/r, 1^, (and (p*) yield the one-particle irreducible vertex 
functions. 



ni. MEAN-FIELD THEORY 

As a warm-up for the renormalization group treatment it 
is instructive to recapitulate the mean-field theory for the su- 
perfluid phase in the functional integral formalism.^'' In mean- 
field approximation bosonic fluctuations are neglected, that is, 
the bosonic field 4> is fixed instead of being integrated over all 
possible configurations. The fermion fields can then be in- 
tegrated exactly. The (fixed) bosonic field is determined by 
minimizing the effective action as a functional of (p. For a 
homogeneous system, the minimizing (pq can be be non-zero 
only for q - Q. We denote the minimum by a. Substituting 
00 ^ Q" + 00 yields 



ipkaHh) - ^k) 4fkcr - a*— a 

ka U 



Jkc 

+ J^i^f-kl^kt a + ^k-\^-ki a") 
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A necessary condition for a minimum of the effective action 
is that its first derivative with respect to (p (or (p*), that is, the 
bosonic 1 -point function r^''(<7), vanishes. In other words, 
terms linear in (or 0*) have to vanish in the effective action. 
For q ^ 0, r^'^(g) vanishes for any choice of c in a homoge- 
neous system. For q - Q and in mean-field approximation, the 
1 -point function is given by 



r<"(0) 



(6) 



where (. . .) denotes expectation values. The first term on the 
right hand side corresponds to the contribution -■iQ'0* to Fq 
in the third line of Eq. while the second term is generated 
by contracting the fermions in the contribution proportional 
to 0* in the forth line of Eq. (|5]). In the absence of bosonic 

fluctuations there is no other contribution to F^''. From the 
conditionF^''(0) = one obtains 



(7) 



which relates c to a fermionic expectation value. 

We now turn to the fermionic 2-point functions. The normal 
fermionic propagator G/o-(^) - -{^'ka-^kcr) and the anomalous 
propagators Ff{k) = -{fpk]^-ki) and F fik) = -{fp-ki^Pk^) can 
be conveniently collected in a Nambu matrix propagator 

The anomalous propagators satisfy the relations Ff{k) - 
F*^{k) and Ff(-k) - Fjik). In (our) case of spin rotation in- 
variance the normal propagator does not depend on cr, and 
therefore, Gjiik) - Gf[(k) - Gfik). 

In mean-field theory, the fermionic 2-point vertex function 
r^^'' = can be read off directly from the bare action in 

the form Eq. (|5]i: 



rf(^) 



ih) - a 
a* iko + 



(9) 



The off-diagonal elements are due to the terms in the second 
line of Eq. (|5]l. Tadpole contributions which are generated 
from the terms in the lines 3-5 of Eq. (|5]) cancel exactly by 
virtue of the condition F^'* = 0. In the absence of bosonic 

(2) 

fluctuations there are no other contributions to Ty . Inverting 

(2) 

T\ and using = ^k from reflection symmetry yields 



Gf(k) 
Ff(k) 



-iko - gk 
k^ + 



k^ + F^ 



(10) 
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where Ek = (^^ + |Ap)'^^ and A - a. We observe that in 
mean-field theory the bosonic order parameter a is equivalent 
to the gap A in the fermionic excitation spectrum. Eq. ^ 
corresponds to the BCS gap equation 



(12) 



We finally compute the bosonic 2-point functions in mean- 
field theory. The bosonic propagators Ghiq) = -{(pq<p*} and 
Fbiq) - -{<Pq<P-q) - -{<P*-q<P*qY form the matrix propagator 



Gbiq) = 



Ghiq) Fh(q) 
Fliq) Gh(-q) 



(13) 



,(2) 



Note that Fbi-q) - Fbiq). The bosonic 2-point function 
is equal to -G^'. We define a bosonic self-energy via the 
Dyson equation (G;,)"' - (Gm))"' - '^b, where the bare prop- 
agator corresponding to the bare action Fq is given by 



Gboiq) 



U 
U 



(14) 



In mean-field theory, only fermionic bubble diagrams con- 
tribute to the bosonic self-energy: 



^^biq) 



K{q) L(q) 
L*{q) Ki-q) 



(15) 
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where 



Gf(k + q) Gf(-k) 



Uq) 



= J^Ff(k + q)Ff(-k) . (17) 



In the absence of bosonic fluctuations, there are no other con- 
tributions to "Lt,. Tadpole diagrams cancel due to F^'^ - 0. 
Note that K(-q) - K*(q) while L(-q) - L(q). Inverting the 
matrix Gho - one obtains the bosonic propagator in mean- 
field approximation 



GM = — 

d(q) 



1 (U-'-K(-q) L(q) 



U{q) - K(q) 



(18) 



with the determinant d(q) ^ \U'^ - Kiq)\^ - \L(q)\^. 

Using the explicit expressions ( fTOl ) and ( fTTT i for G/ and Ff, 
respectively, one can see that 



f/-' - K(0) + |L(0)| 



(19) 



which vanishes if A is non-zero and satisfies the gap equation. 
Hence d{q) has a zero and Gb(q) a pole in q' = 0. This pole 
corresponds to the Goldstone mode associated with the spona- 
neous breaking of the U(l) symmetry of the model. For small 
finite q and qo, the leading ^'-dependences of |f/"' - K(q)\ and 
L(q) are of order |qp and q^. Hence the divergence of Gbiq) 
for — > is quadratic in q and qo. Continuing to real fre- 
quencies one obtains a propagating mode with a linear disper- 
sion relation. As expected, the second pole of d(q) is gapped 
and features a quadratic momentum dispersion. 

By appropriately tailoring the fRG-trunction in Sec. |IV B 21 
to the results of this mean-field calculation in the superfluid 
phase, we incorporate the effects of transversal (Goldstone) 
fluctuations as well as longitudinal (radial) fluctuations into 
our computation. 



IV. FLOW EQUATIONS 

In this section we derive approximate flow equations by 
truncating the exact flow equation^^ for the effective action. 



write down the exact flow equations it is convenient to com- 
bine fermionic and bosonic fields in a superfield S, where 
fermions and bosons are distinguished by a statistics index 
s = b,f, that is, >S/, = O and S/ = ^F. The superpropagator 
G(q) - -{SgSq} is diagonal in the statistics index. 

The flow equation describes the evolution of the effective 
action as a function of a flow parameter A, usually related to 
a cutoff. In this work we use sharp frequency cutoffs which 
exclude bosonic fields with |frequency| < At, and fermionic 
fields with |frequency| < Af from the functional integral. 
Thereby both fermionic and bosonic infrared divergences are 
regularized. Both cutoffs are monotonic functions of the flow 
parameter, Ab{A) and A/(A), which vanish for A ^ and 
tend to infinity for A — > oo. The cutoff can be implemented 
by adding the regulator term 



(?) % + j^'^k Rf(k) *Pi 



to the bare action, where 



R^(^) = [G,o(^)]-' - [Xs(ko)G,oik)]- 



(22) 



(23) 



for s - b,f, andx')^(ko) = ©d^ol - A^). This term replaces the 
bare propagators G^o by G^g = x'^Gso- 

Integrating e"'^""^'^ in the presence of source fields coupling 
linearly to S and S yields the cutoff-dependent generating 
functional for connected Green functions 



-Fo IS.SWR" [S,S]+(S' ,S)+(S,S' ) 



G'^[y,y] = - log I DSDSe 

(24) 

where the bracket (., .) is a shorthand notation for the inner 
product for superfields. The cutoff-dependent effective action 
F^ is defined as 



(25) 



where £.0^ is the Legendre transform of G'^. As a function of 
decreasing cutoff, F'^ interpolates smoothly between the bare 
action Fq for A - oo and the full effective action F recovered 
in the limit A — > 0. 

The evolution of the effective action follows the exact flow 
equation 



A. Exact flow equation 

To simplify our notation, we use fermionic Nambu fields 



dA 



r^[S,S] = Str 



(26) 



and bosonic Nambu fields 



(20) 



(21) 



The fermionic and bosonic matrix propagators are then given 
by Gf{k) = -i^S'k) and Gi(^) = -((D,^<i),^) , respectively. To 



-(2) A 



dSdS 



(27) 



The supertrace Str traces over all indices with a plus sign for 
bosons and a minus sign for fermions. Note that the defini- 
tions of F'^ vary slightly in the literature. In particular, F'^ 
is frequently defined as the Legendre transform of with- 
out subtracting the regulator term 'R'^, which leads to a simple 
additional term in the flow equation. 
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To expand the functional flow equation ( |26] ) in powers of 



the fields, we write the Hessian of T as 



Gf(k) 



Gbiq) 



(28) 



where ^^'^^[5, 5] contains only terms which are at least 
quadratic in the fields. Defining = [(G^) ' - R^] ' and 
expanding in powers of f ^ yields 



-r^ = -Str(R^G^) 

- Str [G;,^(f + f(2)AG^f (2)^ + ■■■)], (29) 



where 



G/A _ C'A-nA 



(30) 



The socalled single-scale propagator G^^ has support only for 
frequencies at the cutoffs, that is, for |A:o| - A,. Expanding 
both sides of Eq. ( |29l ) in powers of the fields and comparing 
coefficients yields a hierarchy of flow equations for the vertex 
functions. 

To describe spontaneous symmetry breaking we expand 
around a generally non-zero value of the ^ = component of 
the bosonic field, that is, we expand r'^[i^, iff, + (f>] in pow- 
ers of i/', if/, (p, (p*, where a'^ may be non-zero. The expansion 
point a'^ will be chosen such that the bosonic 1 -point function 
r^'''^ vanishes for all A. In this way tadpole contributions are 
avoided. The A-dependence of the expansion point a'^ leads 
to terms proportional to a'^ - di^a^ stemming from the total 
A-derivative on the left hand side of the flow equation. 

In the following we will frequently suppress the superscript 
A attached to cutoff'-dependent quantities. 



B. Truncation 

The exact efifective action contains an infinite number of 
terms of arbitrary order in fermionic and/or bosonic fields. 
We now describe which terms are kept and how they are 
parametrized. We keep all terms which are crucial for a qual- 
itatively correct description of the low-energy behavior of the 
system. We distinguish the symmetric regime, where a - Q, 
from the symmetry broken regime, where a Q. The former 
applies to large A, the latter to small A. 



I. Symmetric regime 

The bare action Eq. ^ contains quadratic terms for 
fermions and bosons, and an interaction term where bosons 
couple linearly to a fermion bilinear In the effective action 
we keep these terms with generalized cutoff-dependent pa- 
rameters and add a bosonic self-interaction of order |<^|'*. The 
latter is generated by the flow and becomes crucial when the 
quadratic part of the bosonic potential changes sign. Other 
interactions generated by the flow are neglected. 





FIG. 1 : Diagrammatic constituents of our truncation in tlie symmet- 
ric regime as described in Sec. lIVBTI 



For our choice of parameters (relatively small U), the 
fermionic propagator receives only Fermi liquid renormaliza- 
tions, leading to a slightly reduced quasi particle weight and 
a weakly renormalized dispersion relation. We neglect these 
quantitatively small effects and leave the quadratic fermionic 
term in the action unrenormalized, that is. 



= - I ^kaiih) - ^k) if^ka 
Jko- 



(31) 



corresponding to an unrenormalized fermionic propagator 



Gf(k) = Gfoik) 



1 



(32) 



In the bare action the term quadratic in bosons contains 
only a mass term. In the effective action this mass decreases 
with decreasing cutoff until it vanishes at a critical scale A^, 
which marks the transition to the symmetry-broken regime. 
As the mass decreases, the momentum and frequency de- 
pendence of the bosonic 2-point function becomes important. 
The latter is generated in particular by fermionic fluctuations. 
For small q, the leading q-dependence is of order |qp. The 
leading frequency-dependent contribution to the real part of 
the bosonic 2-point function is of order qfy The frequency- 
dependence of the imaginary part is generally of order ^o, but 
the prefactor is very small, which is related to the fact that it 
vanishes completely in case of particle-hole symmetry. Fur- 
thermore this small imaginary part does not have any qualita- 
tive impact on the quantities we compute in the following. We 
therefore neglect this term and make the ansatz 



where ^ 2 2:ti(l - cos^,) is fixed, while m^, Z/,, and Ah 
are cutoff-dependent numbers. The function a>^ has been cho- 
sen such that the quadratic momentum dependence for small q 
is continued to a periodic function defined on the entire Bril- 
louin zone. The initial conditions for the parameters in the 
bosonic 2-point function. 



Gbiq) 



+ Aba>l + ml 



(34) 



can be read off from the bare action as = \2/U\ and Zj, = 
= 0. 
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The interaction between fermions and bosons remains reg- 
ular and finite near Ac- It can therefore be parametrized as 



and 



^8 I {4f-k+'ii^kVi]4>q + ^k+'i]^-k+iifci) ^ (35) 

JkM 



where the coupling constant g depends on the cutoff, but not 
on momentum and frequency. The initial condition for g is 

^ = 1- 

The flow generates a bosonic self-interaction which plays 
a crucial role near and in the symmetry-broken regime, that 
is, when the bosonic mass term becomes small and finally 
changes sign. The most relevant term is a local |0|"^-interaction 



(36) 



where the coupling constant A depends on the cutoff but not 
on momentum and frequency. The initial condition for A is 
A = 0. 

The propagators and interaction vertices are represented di- 
agrammatically in Fig. 1. 



2. Symmetry broken regime 

For A < Ac the effective action develops a minimum at 
0?=o - a + 0. Due to the U(l) symmetry associated with 
charge conservation the minimum is degenerate with respect 
to the phase of a. In the following we choose a real and 
positive. We decompose the fluctuations of the bosonic field 
around ff in a longitudinal (real) and a transverse (imaginary) 
part, cr and n, respectively: 



where 7r_^ - n* and cr_^ 



(37) 



The bosonic part of the effective action consists of a local 
potential and momentum and frequency dependent contribu- 
tions to the 2-point functions. A U(l)-symmetric local poten- 
tial of order \(pf with a minimum in a has the form 

U^^W - ^ J {m^ - \aff , (38) 

where (p is the original (unshifted) bosonic field. Substituting 
<p — > a + cr + in and expanding in cr and n yields a mass 
term for the cr-field and various interaction terms (see below). 
No mass term for the Tr-field appears, as expected, since the 
TT-field is a Goldstone mode.^- 

The leading momentum and frequency dependence of the 
bosonic 2-point function is quadratic in q and ^o, both for the 
cr- and 7r-component. Hence we make the following ansatz for 
the quadratic bosonic contributions to the effective action 



o-Jm^ + Z, 



(39) 



1 r 2 2 

r;r;r = 2 J ^-^^"1^ + ^JrWg) 



(40) 



where m^-, Z^-, A^-, Z,^, and A„ are cutoff dependent real num- 
bers. The propagators for the cr and n fields thus have the 
form 



1 



OT^ + Za-ql + Aa-Oj\ 



and 



GK(.q) = - 2 \ — 2 ' 



(41) 



(42) 



respectively. The longitudinal mass is determined by the 101"* 
coupling A and the minimum a as 



ml^A\a\^ 



(43) 



The small imaginary contribution to F^.^ generated by 
fermionic fluctuations in the absence of particle-hole sym- 
metry, mentioned already above, gives rise to an off-diagonal 
quadratic term Ya-ji with a contribution linear in ^q. Since this 
term would complicate the subsequent analysis without hav- 
ing any significant effect, we will neglect it. 

The bosonic interaction terms obtained by expanding 
U^'^iip] ai-e 

rCT4 = yo-t I Cr-q-pO--q'+pCrq'Crq , (44) 



(45) 



•Jq,q',p 

- Jn* \ ^-q-p^-q'+pTTq'TT, 
Jq,q',p 

^a-^K^ = Ja-'-n^ I Cr^q^pCr^q^+pTTq^nq , (46) 

Jq,q',p 

r^r-' = 7cr> I 0--q-q'Crq,Crq , (47) 
Jq,q' 

'^crn^ = 7a-,:-- I O-^q^q.Jlq.nq , (48) 

with =y^= Aji, y^2„i = A/4, and y^j = y^^i = Aa/2. 

Let us now discuss terms involving fermions in the case 
of symmetry breaking. In addition to the normal quadratic 
fermionic term F^,;, defined as before, see Eq. dSTT l. the anoma- 
lous term 



{Ai//-kHf'k^ + A*tffk-\tf'-ki) 



(49) 



is generated in the symmetry-broken regime, where |A| is a 
cutoff-dependent energy gap. The phase of A is inherited from 
the phase of a while its modulus is generally different, due to 
fluctuations. Since we have chosen a real and positive, A is 
real and positive, too. The normal and anomalous fermionic 
propagators Gf and Ff corresponding to F^^ and F^^ have 
the standard mean-field form as in Eqs. ( fTUl i and ( fTTI ). with 
- (^2 ^ |A|2)'/2, but now A is not equal to a. 
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Gf(k) = ^ , Ff{k)= » ^ 

GAk) = , GM= 'N.^-v^ 




FIG. 2: Diagrammatic constituents of our truncation in tiie 
symmetry-brolcen regime as specified in Sec. IIV B 21 



C. Approximate flow equations 

Inserting the above ansatz for the truncated effective action 
into the exact flow equation and comparing coefficients yields 
a set of coupled flow equations for the cutoff dependent pa- 
rameters. The various contributions can be conveniently rep- 
resented by Feynman diagrams. The prefactors and signs in 
the flow equations could be extracted from the expansion of 
the exact flow equation, Eq. ( |29] l. However, in practice we 
determine them by comparison to a conventional perturbation 
expansion. 

All contributions to our flow equations correspond to one- 
loop diagrams with only one momentum and frequency inte- 
gration, as dictated by the structure of the exact flow equa- 
tion in the form ( |29l ). One of the propagators in the loop is a 
bosonic or fermionic component of the single-scale propaga- 
tor G^, the others (if any) are components of Gr. 

For a sharp frequency cutoff the frequency variable running 
around the loop is pinned by GJj(A:o) to ko - +A/, or fco = ±^/- 
Hence the frequency integral can be performed analytically. 
The problem that the integrand contains also step functions 
x'^ih)) - ©d^ol - Aj) can be treated by using the identity 
J dx 6{x - X{)) f[x, &(x - xo)] - Jq du f(xo, u), which is valid 
for any smooth function /. 

More specifically, in the present case the one-loop diagrams 
are evaluated for vanishing external frequencies, such that all 
internal propagators carry the same frequency. In loops in- 
volving only either only bosonic or only fermionic propaga- 
tors, one can use the identity 



In addition to the interaction between fermions and bosons 
of the form Eq. (l35T l. an anomalous term of the form 

) (50) 

is generated in the symmetry-broken regime. Inserting the 
decomposition of (p in longitudinal and transverse fields into 
the normal and anomalous interaction terms, we obtain 

r^V = go- (lA-^+axiAt+fTO"? + 'At+fT'A-/t+2iCr-9) > (51) 

Jk,q ^ ^ - . z / 

where gcr = g + g and gn = g - g- Fermions couple with 
different strength to the cr- and 7r-field, respectively. 

A diagrammatic represention of the various progagators 
and interaction vertices in the symmetry-broken regime is 
shown in Fig. 2. We finally note that in a previously reported 
attemptii to truncate the fRG flow in a fermionic superfluid 
the gap A was not taken into account in the fermionic propa- 
gator, and no distinction between longitudinal and transverse 
fields was made for the bosonic Z-factors in the symmetry- 
broken regime. 



n f dkoG'^g(ko)A[GMAr-'^A'^ ^ [G,(^o)A]", 

(53) 

valid for any matrix A, to replace the frequency integration 
by a frequency sum over +As while replacing all the propa- 
gators in the loop by Gj. The factor « corresponds to the n 
possible choices of positioning G'^^ in a loop with n lines, and 
A, - d As Id A. For Ah - Af the above formula holds also for 
the superpropagatorGft, such that is applies also to loops with 
mixed products of bosonic and fermionic propagators. For 
At, + Af, mixed loops contribute only if the single-scale prop- 
agator is associated with the larger cutoff. For example, for 
Af, > Af, the single-scale propagator has to be bosonic, since 
Gb vanishes at |A:o| - Af. On the other hand, for |A;o| - Ah one 
has G fRiko) = G f{ko), and for the integration of the bosonic 
factors in the loop one can again use Eq. ( |53l ). 

For loop integrals involving the frequency sum over +As 
and the momentum integral over the Brillouin zone we use 
the short-hand notation 



r =^ y r 



d^k 
(2^ 



(54) 



where A' - dAs/dA. 
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FIG. 3: Feynman diagrams representing the flow equations in the 
symmetric regime. 



1. Symmetric regime 

Here we choose Kb - Kj - K. One is in principle free to 
choose the fermionic and bosonic cutoff independently. We 
have checked that the concrete choice of A/, and Aj in the 
symmetric regime does not change the final results for A ^ 
much. The Feynman diagrams contributing to the flow in the 
symmetric regime are shown in Fig. [3] 

The flow of the bosonic mass is given by the bosonic self- 
energy at vanishing external momentum and frequency, that 
is, 

Oa'^^S^ [ Gf(k)Gf(-k) + ^ [ GiAq). (55) 

^ Jk\A ^ Jq\h 

The fermionic contribution to ^Am^ is positive, leading to a 
reduction of upon decreasing A, while the bosonic fluctu- 
ation term is negative (since Gbiq) < 0). The flow of Zt is 
obtained from the second frequency derivative of the bosonic 
self-energy as 

d^Zb = ^2 r dlGf(k + q) Gf(-k)\ . (56) 

Jk\A 

Similarly, the flow of is obtained from a second momentum 
derivative of the bosonic self-energy: 

dAAb ^g^ f d\G f{k + q)G /(-fe)l^^o , (57) 

where - 2 '^Ir ^^'^'^^ the bosonic self-energy is 
isotropic in q to leading (quadratic) order in q, the results 
do not depend much on the direction in which the momen- 
tum derivative is taken. The bosonic tadpole diagram in Fig. 3 
contributes only to nib, not to Z;, and Ab, since it yields a mo- 
mentum and frequency independent contribution to the self- 
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FIG. 4: Contributions to the bosonic 1-point vertex and fermion gap 
below Ac. 



energy. Finally, the flow of the \cf)\ coupling is given by 

- -V r [Gf(k)nGf(-k)]^ + ^A^ f [Gb(q)f . 

Jjt|A ^ Jq\K 

(58) 

Within the truncation of the effective action described in 
Sec. IV.B there is no contribution to the flow of the interaction 
between fermions and bosons in the symmetric phase. The 
coupling g remains therefore invariant. 



2. Symmetry broken regime 

In the limit A — > we are forced to choose hj <^ Kb to 
avoid an artificial strong coupling problem, as will become 
clear below. We therefore choose A/ < A;, in the entire 
symmetry broken regime, which implies that the frequencies 
in mixed loops with bosonic and fermionic propagators are 
pinned at the bosonic cutoff. The precise choice of the cutoffs 
will be specified later 

We first derive the flow equation for the minimum of the 
bosonic potential a, which is derived from the condition that 
the bosonic 1-point vertex FJ^^ be zero for all A. The flow 
equation for fJ^'' reads 

d^Y^^^ ^mld^a+2gA Ff(k)+^f [3GAq)+GAq)] ■ 

Jk\Af ^ Jq\\,, 

(59) 

The various contributions are represented diagrammatically in 
Fig. m The first term is due to the cutoff-dependence of the 
expansion point around which the effective action is expanded 
in powers of the fields. The condition cJaF*^* — yields 

d^a = r Ff{k) VhGM) + GAq)] ■ (60) 

Jk\Af 2q' J,^|a„ 

We have used Eq. ( |43] ) to simplify the last term. The fermionic 
contribution to 5aQ' is negative, leading to an increase of a 
upon decreasing A, while the bosonic fluctuation term is pos- 
itive and therefore reduces a. The behavior of Eq. ( I6OI 1 in the 
vicinity of the critical scale, A < Ac, when a and are small 
is shown below in Sec. IV A| 

The flow of A is obtained from the flow of the anomalous 




• -o- 



-i y 



■o- 




FIG. 5: Diagrammatic representation of tlie contributions to the 
bosonic mass in the symmetry-broken regime. 



FIG. 6: Diagrams renormalizing the bosonic Z- and A-factors for 
A < A^. 



component of the fermionic self-energy as 



5aA = d^a 



-j 



Ffiq 



■k)[glGAq)-glGM]\ 



'li:=(0,kf) ■ 

(61) 

The first term, due to the cutoff-dependence of the expansion 
point for the effective action, Hnks the flow of the fermionic 
gap to the flow of the bosonic order parameter The second 
term captures a coiTection to the relation between a and A due 
to bosonic fluctuations as illustrated in Fig.lD 

The flow of the mass of the longitudinal order parameter 
fluctuations (cf. Fig.|5]) is obtained from the self-energy of the 
cr fields at zero momentum and frequency as 



Oa'y = slj^^^ [Gf{k)Gfi-k)-Ffik)Ff(-k)] + 3^d. 
J f [3GAq) + GAq)] 

f [9Gl{q) + Gliq)] . (f 



(Aaf 



The second term in this equation is due to a product of the 3- 
point vertex y^-i and (9a o- arising from the cutoff dependence 
of the expansion point for the effective action. The flow of A 
can be computed from the flow of and a via the relation 
Eq. A = ml/\a\^. 

The flow of Zcr is obtained from the second frequency 
derivative of the self-energy of the cr fields, which yields 



5aZ^ ^ glf dl[Gf(k + q)Gf(-k)-Ff(k + q)Ff(-k) 
Jk\Af 



\q=0 



f dl [9GJk + q) G^k) + GJk + q) G„{k)] I . 

(63) 



The flow of Acr is given by the same equation with 5^^^ replaced 
by 4 



In the flow of Z„ there are strong cancellations of differ- 
ent terms originating from bosonic fluctuations. These can- 
cellations are related to Ward idenities which guarantee that 

remains finite such that the Goldstone mode is not renor- 
malized substantially. Hence, we keep only the fermionic 
fluctuations, that is, 

5aZ, ^gl \ d^iGfik + q)Gf(-k) + Ff(k + q) Ff(-k)]\ . 

Jk\Af 

(64) 

For the flow of A„ we obtain the same equation with d^^^ re- 
placed by 5q. The terms contributing to the flow of the Z- and 
A-factors are illustrated diagrammatically in Fig.|6] 

In the symmetry broken regime, there are also contributions 
to the flow of the interaction between fermions and bosons due 
to vertex coiTections with bosonic fluctuations (see Fig. |2l), 
yielding 



dAg^ ^8^ f \Fj(k -q)- \Gf{k - q)f] 



*=(0,kf) 



and 



/ 

Jq\Ai 



[glGM)-glGM)\ , 



F]{k-q) + \Gf{k-qf 



[glGAq)-glGM)\ 



(65) 



A=(0,kf) 



(66) 



The right hand sides are dominated by the contribution from 
the n propagator, which tends to reduce go-, gn for decreasing 
A. 



D. Relation to mean-field theory 

Before solving the flow equations derived above, we first 
analyze what happens when contributions due to bosonic flue- 



10 







FIG. 7: Fermion-boson vertex corrections below A,. 



tuations are neglected, and relate the reduced set of equations 
to the usual mean-field theory (Sec. 111). 

In the absence of bosonic fluctuations, g - go- - g„ = I- 
Furthermore, the bosonic order parameter a and the fermionic 
gap A are identical: a = A. Since the bosonic cutoff is irrele- 
vant here, we can choose A/ = A. 

The flow equation for the bosonic mass in the symmetric 
regime, Eq. ( |55] l, simplifies to 

m? f 

5a ^= Gf(k)Gfi-k), (67) 

^ Jk\A 

where Gf{k) - Gfo(k) = (iko - ^k) This equation can be 
easily integrated, yielding 



It is complemented by the flow equation for the mass of the cr 
field, Eq. ( |62] |. which becomes 

5a^ = j^^^[\Gf(kf -Fj(k)] + 3y^.dAA (70) 

in the absence of bosonic fluctuations, with ja-^ = AA/2 = 
OT^/(2A). The propagators Gj and Fj have the usual BCS 
form, as in Eqs. ( fTOl i and (fTTT i. 

A numerical solution of the coupled flow equations ( |69l ) 
and dTOl l yields a gap A which is a bit smaller than the BCS 
result obtained from the gap equation ( fT2b . The reason for 
this discrepancy is the relatively simple quartic ansatz ( l38l l for 
the bosonic potential. The complete bosonic potential is non- 
polynomial in \<p\^ even in mean-field theory. Restricted to 
the zero momentum and frequency component of <p it has the 
formal 



\U\ 



-I- 



(71) 



The kernel of the 3 -point vertex ja-i obtained from an expan- 
sion of this mean-field potential around a finite order parame- 
ter A reads 

= -2 J [\F]{k) - Ff(k) \Gf(k)f-] (72) 

at zero frequencies and momenta. Inserting this into (iTOb . the 
flow of can be written as a total derivative 

dA^^-d^f f\\Gf(k)\^-Fj(k)\, (73) 

^ J|*ol>A -Jk 

where the A-derivative on the right hand side acts also on A, 
generating the term proportional to y^^i ■ Integrating this equa- 
tion with the initial condition m^r = at A = Ac , yields 



(74) 



which is the correct mean-field result. With OTo- given by ( 1741 ), 
the flow equation ( |69t yields the correct mean-field gap. The 
easiest way to see this, is to write the BCS gap equation in the 
presence of a cutoff in the form 1 = J|j. |>aX^ ' ^/(^)' 
and take a derivative with respect to A. 

It is instructive to relate the above flow equations for A and 
nicr to the flow equations for the BCS mean-field model ob- 
tained in a purely fermionic RG.^ For a sharp frequency cut- 
off, those flow equations have the form 



5aA = -(V+W) 



m. 



L ^^^^^ ' 



(75) 



\U\ XoI>aX^^+^^' ^^^^ dj^(V+W) ^ (V+Wfd^[ [\\Gf(k)f-Fj(k)],a6) 

J\ko\>A Jk 



nil, vanishes at a critical scale A^ > 0. 

The flow equation for A (= a) in the symmetry broken 
regime A < A^ , Eq. i6T[ . is reduced to 



vA=-^ r . 

mfr Jk\A 



Ffik) 



(69) 



where V is a normal two-fermion interaction in the Cooper 
channel, while W is an anomalous interaction corresponding 
to annihilation (or creation) of four particles. With the identi- 
fication -V + W these equations are obviously equiv- 
alent to f6% and ( |73] |. The above flow equation for V + W 
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is obtained from a one-loop truncation complemented by ad- 
ditional self-energy insertions drawn from higher order dia- 
grams with tadpoles. ■'■^ These additional terms correspond to 
the contractions with 7,^.3 in the present bosonized RG. 



V. RESULTS 



--0- 



FIG. 8: Goldstone fluctuations determining the infrared asymptotics 
A. Flow for A < A, for A 0. 



For A slightly below Ac, the flow equations can be ex- 
panded in the order parameter To leading order, the order 
parameter a and the gap A are identical, ff = A . The fluctua- 
tion term in the flow equation ( I6TI 1 for A is quadratic in a, and 
also the flow of g^r yields only corrections beyond linear order 
to the relation between a and A. 

Near Ac, the flow equation (I6OI 1 for a can be written as 



[3GAq) + GAq)], 



(77) 



Note that we 



where I - J^|^ (^q +^^) ' , evaluated for A = A^ 

have replaced the ratio rn^la^ by A in the first term on the right 
hand side of the flow equation. Integrating the flow equation 
one obtains 



-1+ { ['iGAq) + GM'\ 



{Ac - A) (78) 



A=A, 



for A < Ac. The order parameter a and the fermionic gap are 
thus proportional to (Ac - A)'^^ for A < Ac. 

Inserting Eq. ( |78] l into the flow equation (|62| i for and 
neglecting the last fluctuation term, which is of higher order 
in a, one obtains 



dAm, 



Jq\Kh 



This shows that the flow of a and is indeed consistent 
with the relation = Aa^ following from the ansatz for the 
bosonic potential. 



B. Infrared asymptotics 

In the infrared limit (A 0) the key properties of the flow 
can be extracted from the flow equations analytically. The be- 
havior of the bosonic sector depends strongly on the dimen- 
sionality of the system. We consider dimensions d >2, focus- 
ing in particular on the two- and three-dimensional case. 

The bosonic order parameter and the fermionic gap saturate 
at finite values in the limit A — > 0. The fluctuation correc- 
tions to Sao; and cJaA involve the singular Goldstone propaga- 
tor G„ only linearly and are therefore integrable ind > 1 . The 
fermion-boson interactions g^r and g„ also saturate. The finite- 
ness of Z„ and A„ is guaranteed by Ward identities.-^ In our 
truncation Z„ and A„ remain finite since only fermionic fluc- 
tuations contribute to the Goldstone propagator. We choose 



Afe = A in the following. The choice of Af (as a function of 
A) will be discussed and specified below. 

The flows of m^. A, Za-, and A,j- are dominated by terms 
quadratic in G„ for A — > 0, see Fig. [8] Using - Aa^, we 
obtain the asymptotic flow equation for A from Eq. ( |62] | in the 
simple form 



^^A = A- 



j' 

Jq\h 



Gl(q) . 



(80) 



The above integral over is proportional to A"'"'* for small A 
in dimensions d < 4, implying that A scales to zero in li < 3. 
In two dimensions one obtains Gliq) - ^ A~^ for 

A/ A obeys the 



'q\A 

small A, such that the rescaled variable A 
flow equation 



dA 



t/log A 



-A + 



(81) 



which has a stable fixed point at A* 
bosonic self-interaction vanishes as 



'47:^A„Z„ ' 

An^A^Zn. Hence, the 



for A ^ 



(82) 



in two dimensions. Consequently, also the radial mass m^. of 
the Bose fields vanishes linearly in A. In three dimensions 
(79) one has J^^^ G^{q) oc A"' for small A such that A and scale 

to zero logarithmically for A — » 0. Since is the domi- 
nant contribution to the denominator of Go- at small momenta 
and frequencies, the scaling of to zero as a function of A 
imphes that Go- (at A = 0) diverges as 



Ga-(sq) oc 



for d — 2 



Ga-(sq) oc log s foTd-3 



(83) 
(84) 



in the limit 5 — > 0. Although derived from an approximate 
truncation of the functional flow equation, this result is exact 
even in two dimensions, where the renormalization of is 
very strong. This is due to the fact that the scaling dimen- 
sion of is fully determined by the scaling dimension of the 
Goldstone propagator and the existence of a fixed point for A, 
but does not depend on the position of the fixed point.^^ 
The flow of Zo- is given by 

5aZ^ = ^ dlGAp + q)GAq)\,,_, (85) 

^ Jq\A 

for small A. The integral over the second derivative of G;rG;r 
is of order A''"^. In two dimensions the coupling A vanishes 
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linearly in A, such that 5aZ(^ oc A ^, implying that Zg- di- 
verges as A"' for A — > 0. Hence, the term Z^-^q with Iq^l - A 
in the denominator of Go- scales linearly in A, as m^. In three 
dimensions f^^^d^pfijiip + q)G„(q)\p=o diverges as A"^, while 

A vanishes only logarithmically. Hence 5aZo- oc (log A)"^A"^, 
which is larger than in two dimensions. Integrating over A 
one finds Z,r oc (Alog A)"^, which means that Zo-q^ vanishes 
as (log A)"^ in the infrared limit. This yields a subleading log- 
arithmic correction to the mass term in the denominator of 
Go-. An analogous analysis with the same results as just ob- 
tained for Z,j- also holds for the momentum renormalization 
factor Aa-. A strong renormalization of longitudinal correla- 
tion functions due to Goldstone fluctuations appears in various 
physical contexts. Recently, a singular effect of Goldstone 
fluctuations on the fermionic excitations in a superfluid was 
found in Gaussian approximation.-^^ This singularity appears 
only after analytic continuation to real frequencies, and its fate 
beyond Gaussian approximation remains to be clarified. 

Since m^^ and A scale to zero in the infrared limit in li < 3, 
all purely bosonic contributions to the eff'ective action scale to 
zero. On the other hand, the fermion-boson coupling remains 
finite. One is thus running into a strong coupling problem, in- 
dicating a failure of our truncation, if fermionic fields are inte- 
grated too slowly, compared to the bosons. The problem man- 
ifests itself particularly strikingly in the flow equation for the 
order parameter, Eq. ( |60] |. in two dimensions. Since oc A/, 
for small A^, the fermionic contribution to ^Aff is of order A"' 
if one chooses A/ = A/, = A, leading to a spurious divergence 
of a for A — > 0. 

The problem can be easily avoided by integrating the 
fermions fast enough, choosing A/ «; A^ in the infrared limit. 
In our numerical solution of the flow equations in the follow- 
ing section we will choose Ah - A and A/ = A^/A^ for 
A < Ac, which matches continuously with the equal choice 
of cutoffs for A > Ac. The fermionic contribution to Oao: in 
Eq. ( |60l l is then finite for small A, since the factor A'j = 2A in 

J^l^ compensates the divergence of m^^ in front of the inte- 
gral. Since the fermions are gapped below A^., one could also 
integrate them completely (set A/ to zero), and then compute 
the flow driven by A;, only. The freedom to choose fermionic 
and bosonic cutoffs independently was exploited also in a re- 
cent fRG-based computation of the fermion-dimer scattering 
amplitude in vacuum.-'^ 



C. Numerical results in two dimensions 

In this section, we present a numerical solution of our flow 
equations from Sec. IIVCI in two dimensions. Technically, 
we employ a fifth-order Runge-Kutta integration routine^^ to 
solve coupled, ordinary differential equations. At each incre- 
ment of the Runge-Kutta routine, two-dimensional integra- 
tions over the whole Brillouin zone have to be executed. For 
this purpose, we employ an integrator for singular functions^^ 
with relative error of less then 1 %. In particular, for A ^ A^ it 
is imperative to operate with sufficiently accurate routines as 
the integrands are large and small deviations result in a signif- 
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FIG. 9: (color online) (a): Flows of the bosonic mass, for A > 
and for A < A^. (b): Quartic coupling A. 



icant spread in the final values for A ^ 0. 

We fix our energy units by setting the hopping amplitude 
t - I. We choose a chemical potential yU = - 1 .44 correspond- 
ing to an average electron density of 1 /2 (quarter-filled band). 
This choice represents the generic case of a convex Fermi sur- 
face remote from van Hove singularities. The only varying 
input parameter is the Hubbard U, which determines the ini- 
tial value of the bosonic mass via = |2/t/|. 

Initially, the flow starts in the symmetric regime with A = 
Ao = 100, where Eqs. (|55]-|58]) determine the evolution. The 
critical scale is determined by the condition m^XAc) - 0. In 
the symmetry -broken regime (A < A^), Eqs. (|60l-l66l) deter- 
mine the evolution. 

In Fig. |9] characteristic flows of the bosonic mass and the 
quartic coupling are shown in double-logarithmic plots for dif- 
ferent choices of the Hubbard U. The sharp de- and increase of 
the bosonic mass marks the region around A^. For small A, the 
flow reaches the infrared asymptotic regime (see Sec. IV Bl) . 
The scale Ajr at which this scaling sets in decreases for de- 
creasing U. The numerically obtained scaling m^, /} oc A is 
consistent with the analytical result of Sec. IV B| 

In Fig.[TO](a) we compare the fermion gap with the bosonic 
order parameter (final values at A = 0) and the critical scale as 
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(b) 
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FIG. 10: (color online) (a): Fermion gap A, order parameter a, and 
the critical scale versus Hubbard U. (b): Exemplary flows for A 
(triangles) and a (lines) each corresponding to one point in (a). 
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FIG. 11: (color online) Flows of Zcr (lines) and Z„ (crosses). 




FIG. 12: (color online) Flows of fermion-boson vertices, g^r (lines), 
and g„ (lines with crosses), for A < A^.. 



a function of U. In Fig.[TO](b) the flow of A and a as a function 
of A is shown for various choices of U. We observe A, a cc 
(Ac - A)'^^ for A < Ac as derived below Eq. (ITST l. The ratio 
A/ Ac, where A is the final gap for A ^ 0, is approximately 1 .4 
for the values of U studied here. As a result of fluctuations, 
the gap is reduced considerably compared to the mean-field 
result 

« 0.25 (86) 

Abcs 

for 2 < U < 4. The main reduction here stems from the 
bosonic self-interactions in the symmetric regime leading to 
a substantial decrease of Ac via the second term of Eq. f55[ . 
A reduction of the gap compared to the mean-field value is 
generally present even in the weak coupling limit f/ — > 0. In 
fermionic perturbation theory second order corrections reduce 
the prefactor of the BCS gap formula even for f/ — » 0.^^ The 
reduction obtained here is slightly stronger than what is ex- 
pected from a fermionic renormalization group calculation.^ 

For A < Ac, Goldstone fluctuations slightly reduce a via 
the term involving G„(q) in Eq. (|60] l. On the other hand, 
the term due to the Goldstone mode in Eq. (|6T| i enhances the 
fermionic gap A relative to a, such that A is generally slightly 



larger than a. This difference will become larger upon in- 
creasing the interaction strength as one enters the regime of 
a Bose gas made of tightly bound fermions. Here, however, 
for relatively weak interactions, the impact of Goldstone fluc- 
tuations on both a and A is very modest. By contrast, the 
impact of Goldstone fluctuations is known to be dramatic at 
finite temperatures (not treated here) in two dimensions, since 
they drive the order parameter to zero.^^ 

In Fig.im we show flows of the Z-factors of the cr- and 
TT-field, respectively. In the symmetric regime, the evolution is 
independent of U. At Ac, the 0-field splits into the cr- and n- 
modes with Zo- diverging in the limit A — > as Z,r A"' (cf. 
Sec. IV Bl ). The Z-factor of the Goldstone field saturates for 
A <^ A. The flows for the A-factors (not shown) parametriz- 
ing the momentum dependence of the cr- and .Tr-propagators 
exhibit very similar behavior. 

Finally, in Fig.[T2]we show flows of the fermion-boson ver- 
tices gu- and g„ for A < Ac. Their relative changes are only of 
the order of a few percent (note the scale of the vertical axis) 
with go- being a bit larger than g„. 
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VI. CONCLUSION 

Truncating the exact fRG flow, we have derived approxi- 
mate flow equations which capture the non-trivial order pa- 
rameter fluctuations in the superfluid ground state of the at- 
tractive Hubbard model, which has been chosen as a pro- 
totype model for attractively interacting fermions. The su- 
perfluid order parameter is associated with a bosonic field 
which is introduced via a Hubbard-Stratonovich decoupling 
of the fermionic interaction. Below a critical scale Ac, the 
bosonic effective potential assumes a mexican hat shape lead- 
ing to spontaneous symmetry breaking and a Goldstone mode. 
The bosonic order parameter is linked to but not equivalent to 
a fermionic gap. The fermionic gap is significantly smaller 
than the mean-field gap, mostly due to fluctuations above the 
scale Ac. Transverse order parameter fluctuations (Goldstone 
mode) below Ac lead to a strong renormalization of radial 
fluctuations. The radial mass and the bosonic self-interaction 
vanish linearly as a function of the scale in two dimensions, 
and logarithmically in three dimensions, in agreement with 
the exact behavior of an interacting Bose gas.^^ On the other 
hand, the average order parameter, the fermionic gap, and 
the interaction between fermions and bosons are aff'ected only 



very weakly by the Goldstone mode. 

Supplementing the flow equations derived above by a shift 
of the chemical potential, to keep the density fixed, one may 
also try to deal with larger values of U . Eagles^^ and Leggett— 
have shown that already the BCS mean-field theory captures 
many features of the condensed Bose gas ground state made 
from strongly bound fermion pairs in the limit of strong attrac- 
tion. Beyond mean-field theory, the difference between the 
fermionic gap and the order parameter a increases at larger 
U. It will also be interesting to extend the present analy- 
sis to r > 0, in particular in view of the possibility of a 
finite fermionic gap in the absence of long-range order in a 
Kosterlitz-Thouless phase at low finite temperatures. 
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